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Abstract 



The entropy current is constructed phenomenologically order by order in the derivative expansion 

On 

^1 | dence, the entropy current of the fluid is mapped to a vector density associated with the event 

■ horizon of the dual geometry. In this work we consider the local horizon entropy current for higher- 



(N 



X 



In the hydrodynamic regime of field theories the entropy is upgraded to a local entropy current, 
re entropy current is constructed phenomenologically order by order in the derivative expansion 
by requiring that its divergence is non-negative. In the framework of the fluid/gravity correspon- 



curvature gravitational theories proposed in arXiv: 1202.2469, whose flux for stationary solutions 
is the Wald entropy. In non-stationary cases this definition contains ambiguities, associated with 
absence of a preferred timelike Killing vector. We argue that these ambiguities can be eliminated 
in general by choosing the vector that generates the subset of diffeomorphisms preserving a nat- 
ural gauge condition on the bulk metric. We study a dynamical, perturbed Rindler horizon in 
Einstein-Gauss-Bonnet gravity setting and compute the bulk dual solution to second order in fluid 
gradients. We show that the corresponding unambiguous entropy current at second order has a 
manifestly non-negative divergence. 

PACS numbers: 04.70.-s, ll.25.Tq, 47.10.ad 
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I. INTRODUCTION 



According to the holographic principle [l, 2|, quantum gravitational theories are equiv- 
alent to certain non-gravitational field theories living in one lower spatial dimension and 
defined on a boundary surface in the higher dimensional bulk spacetime. The only concrete 
realization of holography that we currently possess is based on the AdS / CFT correspondence 
between conformal field theories and their various deformations and quantum gravity (string 
theory) on gravitational backgrounds with negative cosmological constant (for a review see 
[2!). The mysterious nature of holography in general is a crucial aspect of the puzzle of quan- 
tum gravity. On a less fundamental level, holography and the AdS/CFT correspondence 
also offer new ways of investigating unresolved issues both in field theory and in gravitation. 

One important example of this is the hydxodynamic regl me of held theo ry fl, and the 
fluid-gravity correspondence |5|], which originally followed as a special case of AdS/CFT. The 
hydrodynamics associated with a thermal state in the gauge theory is equivalent to the long 
wavelength, long time dynamics of black hole (brane) solutions in the bulk gravity theory. 
One can explicitly construct perturbed black hole solutions order by order as an expansion 
in derivatives of the fluid velocity and temperature and find that the subset of Einstein 
equations constraining data on the boundary surface are the Navier-Stokes equations. The 
essential ingredients needed to relate fluids to gravity are the existence of a horizon in the 
gravitational background that is related to a thermal equilibrium state in the field theory, 
and a derivative expansion around it jfi-sj]. Thus, one can define the relation between fluids 
and gravity on more general backgrounds, for instance the Rindler geometry [9 12] . 

On the gravitational side of the duality, the correspondence has motivated new studies 
of black hole entropy in a dynamical setting, where the horizon surface evolves in time. 
In hydrodynamics the relevant quantity is a local entropy current. In a regime sufficiently 
close to equilibrium, the fluid-gravity mapping implies that the entropy current of the fluid 
low can be constructed in terms of the event horizon geometry as the Bekenstein-Hawking 

in nn 

13. 14] area current |15|, [16(. The thermodynamical Second Law enforcing the positivity 
of the entropy current's divergence is equivalent to Hawking's area theorem in classical 
General Relativity. Further away from equilibrium it has been suggested that the correct 
hypersurfaces in the bulk on which to build the entropy current should be the quasi-locally 
defined apparent horizons (see for example, 17H19|). which also obey the Hawking theorem. 
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The main goal of this paper is to explore within the fluid setting a related issue in semi- 
classical gravity, which is the nature of dynamical horizon entropy in a higher curvature 
theory of gravity. In General Relativity, black hole thermodynamics allows one to clearly 
identify the entropy associated with equilibrium processes as proportional to the cross- 
sectional area of the event horizon. Outside the equilibrium setting, the key requirement for 
defining an entropy current phenomenologically is that it is consistent with the (generalized) 
Second Law. Thus, the link between entropy and area still seems to be robust due to 
Hawking's area theorem the ambiguity is which horizon surface in the bulk is the appropriate 
holographic surface in general. 



In higher curvature theories of gravity the situation becomes more complicated. Wald [20 1 
studied quasi-stationary processes in a general diffeomorphism invariant theory of gravity 
and was able to derive a general formula for the entropy. The relevant quantity is the 
antisymmetric Noether potential Q AB associated with diffeomorphisms along a vector £ A . 
We will focus on the case where the gravitational Lagrangian depends onl y a lgebraically on 



the Riemann tensor £(gAB, Rabcd)- Here the potential has the form 21 



22 



Q AB = J=g (-2C ABCD V C £ D + A£ D V C C ABCD ) , (1) 



with C ABCD = dC/ BRabcd 34J- When the black hole solution is stationary, £ A is the 
timelike Killing vector. The total entropy of the horizon is unambiguous and is proportional 
to an integral of Q AB over any horizon cross-section 



Swald — -j, J Q AB dY>AB ■ (2) 

The bifurcation surface where £ A vanishes and V\a£b] = £ab is a convenient choice for 
actual computations 35] . In Einstein gravity the Wald formula reduces to the Bekenstein- 
Hawking area entropy, but in general the entropy will depend on the both the intrinsic and 
the extrinsic geometries of the horizon surface. 

When the horizon is dynamical, the Wald formula is subject to ambiguities because there 
is no longer a preferred choice of the Killing vector and because one is free to add total 



divergence to the Lagrangian, symplectic potential, and (CD) itself 2l|, |23[. In the paper 
we address the ambiguity in the definition of the vector £ A , and assume that the Noether 
potential has the standard form (CQ) used by Wald. 

Furthermore, there is no known analog of the area theorem in a general higher curvature 
theory, so the Second Law provides no guidance on how to appropriately define an entropy. 



In the context of the fluid- grav ity correspondence, a definition for a local entropy current 



has been recently proposed 



24J , which in a coordinate free form reads 



S A = ^L Q AB £b (3) 
K 



This is to be evaluated on the horizon, so £b is the normal to the horizon and k is the 
surface gravity associated with this normal vector, defined in general as a measurement of 
its non-affinity, I b VbI a = n£ A . Thus the current is effectively a flux of the Noether potential 
through the horizon surface. This expression contains an ambiguity in the definition of the 
vector £ A off the horizon surface in the bulk. In particular the current s A depends on the 



derivative of £ A with respect to the bulk radial coordinate. In 24| it was shown that this 
ambiguity could be eliminated in Einstein's gravity by requiring the vector field to satisfy a 
"weak Killing condition" at the horizon 

£ A (V A £ B + V B £ A ) = • (4) 

When this condition is imposed, the current in ([3]) reduces to the Bekenstein-Hawking 
area current. In higher curvature theories, there was apparently no general way to resolve 
the ambiguity. However, working to first order in fluid gradients, it was shown that the 



weak Killing condition (imposed at this order) again leads to a non-ambiguous current [24 1. 
This current was constructed explicitly in the case of a charged black brane background 
in Einstein-Maxwell theory with a mixed gauge-gravitational Chern-Simons term and its 



divergence was shown to be non- negative 24 1. 



An important question is, whether the ambiguities in ([3]) can continue to be eliminated 
at higher orders in fluid gradients, where the theory is further away from the equilibrium 
state. We propose that this can be naturally done by imposing a more general condition on 
£ (of which the weak Killing condition (j4j) is just one). This is 

£l9rA = , (5) 

where the coordinates are X A = (r, x M ) and £i is the Lie derivative in the direction of £. 
x M are the gauge theory coordinates and r is the holographic radial coordinate. This choice 
follows from the standard gauge imposed on the bulk metric in the fluid-gravity correspon- 
dence g rr = 0, g r ^ ~ u^, where w M corresponds to the fluid velocity. The construction of 
the bulk metric in a derivative expansion is carried out such that this gauge condition is 



preserved. With this gauge choice, one has a map between the boundary and the horizon 
coordinates, and the hydrodynamics is realized in the space. The conditions (jSJ) that 
we propose determines £ A , such that it is preserves the gauge condition. This is, in fact, 
a necessary condition if we want the entropy current constructed from the horizon data to 
be compatible with the hydrodynamics determined by the boundary stress-energy tensor. 
However, while this choice removes the ambiguity in the definition of £ A , we do not know 
whether it guarantees that the entropy current has a non-negative divergence. 

In this work we will investigate the nature of the entropy current in a particular setup. 
We will first construct the full second order solution to Einstein- Gauss-Bonnet gravity (in 
spacetime dimension greater than four) using the recently discovered fluid/Rindler corre- 
spondence 9Hl2| . This correspondence is based on the fact that the Rindler wedge of flat 
Minkowski space acts as finite temperature thermal state analogous to the black brane 
in AdS. Furthermore, like the black brane, the Rindler acceleration horizon has a planar 
topology which allows as discussed above, to make the same controlled expansion in fluid 
gradients as in the fluid-gravity correspondence. Using this machinery, relativistic solutions 



to the vacuum Einstein equations have been constructed to second order [Hi, Il2fl . These 
bulk solutions are dual to a fluid system living on an arbitrary timelike surface S c of fixed 
radial coordinate r = r c . 

So far this is a holographic mapping between two classical theories, but it hints towards a 
full duality between some quantum field theory on S c and the interior region of the Rindler 
geometry. While the nature of holography in asymptotically flat spacetimes is a mystery, it 
is possible that some information can be gleaned from the non-standard properties of the 
dual fluid. In particular, the dual fluid thermodynamics is characterized by zero equilib- 
rium energy density even though there is non-zero temperature. The corresponding viscous 
hydrodynamics of the system is perfectly well-defined, but has the intriguing property that 
only transport coefficients at second (and higher) order in the derivative expansion are af- 



fected by higher curvature terms in the gravitationa 



theory. This implies the shear viscosity 



to entropy density ratio of the fluid is universal 25[. Exploring the hydrodynamics of this 
fluid in more detail is of interest. 



Following the procedure developed in [ll|, [12J, we compute the vacuum solution to second 



order. Since the Rindler background is flat, some simplifications occur and we are able to 
obtain the Gauss-Bonnet corrections to the metric and the fluid stress tensor. This extends 
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the earlier results of 



25| in the non-relativistic limit to the fully relativistic fluid case. In 



particular, we will find to first order in fluid gradients the extension of the horizon normal 
into the bulk, 

and show that the entropy divergence at second order is non-negative. 

The plan of this paper is as follows. In Section II we will review the correspondence 
between hydrodynamics and gravity and motivate the entropy current formula ([3]) in more 
detail. We discuss the nature of its possible ambiguities and why (jSJ) is a natural condition to 
fix I a- In Section III we construct the solution to vacuum Einstein-Gauss-Bonnet gravity up 
to second order in fluid gradients by perturbing around the Rindler background and find the 
holographic fluid stress tensor at second order. In Section IV we compute the entropy current 
using ((3j) and show that its divergence is non-negative. In the discussion, we examine the 
implications of our results and possible extensions to a wider class of examples in the fluid- 
gravity paradigm and beyond. Finally, the appendices contain a more detailed discussion of 
ambiguities and specific calculations of the entropy current. 



II. THE ENTROPY CURRENT 



We begin with a brief review of how the hydrodynamics of a fluid system can be en- 
coded in a gravitational solution in one higher dimension. The key underlying concept in 
hydrodynamics is the notion of local thermodynamic equilibrium. The fluid is described 
by a finite set of macroscopic parameters, which are functions of space and time that vary 
slowly throughout the system, so that in the neighborhood of each point there is an approx- 
imate notion of thermodynamic equilibrium. Thus, a relativistic fluid is characterized by a 
four-velocity and the thermodynamic variables, energy density p, pressure p, tempera- 
ture T, and entropy density s. If the fluid also possesses additional conserved charges, one 
includes the corresponding charge densities and chemical potentials in its description. The 
various thermodynamic variables are related by the equation of state and by the standard 
equilibrium thermodynamical identities. 

Hydrodynamics is an effective description valid at scales L 3> £ m f p , where t m f p is the mean 
free path of the system, determined by the temperature and the nature of the microscopic 
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(field) theory. This means the gradients of all the fluid parameters must be small compared 
to the scale set by the mean free path and thus hydrodynamics is characterized by an 
expansion in derivatives of the fluid variables. At zeroth order, the fluid is in equilibrium 
and the entropy current su^ is conserved. Higher orders in derivatives correspond to viscous 
corrections, which are in general associated with dissipation and increase the entropy of the 
system. 

On the gravity side, we consider a {d + 2) dimensional spacetime and denote the bulk 
coordinates by X A , where the index A runs from (0..d + 1). Typically, one decomposes 
X A = (r, where r is the holographic radial coordinate and x M are the coordinates 
in the field theory/fluid. The index /i runs from (0...d) and therefore d is the number 
spatial dimensions of the fluid system. The metric ansatz corresponding to a fluid in global 
equilibrium is 

ds 2 = g AB dX A dX B = k(r)u^u u dx»dx u - 2u l ,dx tM dr + f{r)P^dx^dx v . (7) 

Here, u^u^ = — 1 and can be thought of as the velocity of a boost. P^ w = h^ u + u il u v 
is the projector orthogonal to and the metric (possibly curved) on which the fluid 
system lives. The functions f(r) and k(r) are determined by the field equations. There 
is an event horizon located at radius such that k(r = r^) = 0. At this location, the 
Eddington-Finkelstein like coordinates chosen for (J2J) are regular. To see that this metric 
corresponds to a holographic fluid, one can compute the Brown- York stress tensor for a 



surface of constant r and show it has the form of a perfect fluid 



26] 



with Hawking temperature 



Tpv = pu^u u + pP^ , (8) 



Ml (q) 

r=r h ■ W) 



4lT 

Note, that in order to describe charged fluids, extra gauge fields are needed in the ansatz. 

The extension to an arbitrary fluid state is straightforward: one promotes m m (x m ), k(r, x M ) 
and f(r, x^ 1 ). The metic ansatz is no longer an exact solution to the field equations, but one 
can work order by order in an expansion in derivatives of these variables as done in j^J. The 
details of this construction for the Rindler metric will be described in the next section. Here 
we note that the perturbed metric solution implies the event horizon is dynamical and its 
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location varies in time and space, r/ l (x M ). The horizon location is determined by solving the 
equation for a null hypersurface 

g AB d A (r - r h {x»)d B {r - r h (z")) = , (10) 

order by order in the derivative expansion. The horizon normal vector £ A then follows from 
£ A = g AB d B {r -r h (x")). 

When working with the the entropy current ([3]) it will be useful consider a coordinate 
gauge adapted to the horizon X A = (f, x M ), where the horizon is always located at zero 
radius f = 0, i.e. 



r = r — Th{x) . (11) 



In these coordinates 



f = g f %=o = 

= f 1f=O > (12) 

and the entropy current ([3]) reduces to 

s A = (0,sn, (13) 

where 

2n 2n - 

Kj Kj 

In the case of Einstein's gravity, where £ = y/—gR (we use units where IQ-kG = 1), 

C ABCD = \{g Ac g BD -g AD g BC ) , (15) 

and the second covariant derivative term in (fill) vanishes identically. In equilibrium one can 
use the Killing equation V a^b = —^b^a to set Vf£ M = — V^f. As a result, ( TH|) reduces to 
the Bekenstein-Hawking area current 

8& R = , (16) 

which for the metric ansatz above (|7|) reduces to s M = 47t/(0) 3/ ' 2 m ai . 

In the dynamical case, there is no longer a Killing vector, but it is possible to use the 
freedom in the radial derivative of Ia (now thought of as a generalization of the Killing 



vector in the bulk which becomes the null normal when evaluated on the horizon) to impose 



the "weak Killing condition' 1 



24] 



\V A £ B + V b £a) = . (17) 



Since the \i component of the equation turns out to be an identity, this amounts to one 
condition 

£ B (V f £ B + Vs4) = \v f {£ B £ B ) + = , (18) 

When evaluated on the horizon this condition relates the surface gravity at higher orders to 
the radial derivative of the norm of I . With this condition imposed, £ M V f £ M = — i^V^if , 
we get that definition (fT4j) again reduces to the Bekenstein-Hawking current. Note, that at 
this stage the derivative expansion has not played a role. As long as weak Killing condition 
can be enforced in a general dynamical situation, (J3J) always reduces to the area current. 

In a higher curvature theory, one must restrict to the fluid-gravity setting and analyze 
the nature of the ambiguities order by order in the derivative expansion. This was done in 



24j for a generic higher curvature theory to first order in fluid derivatives. The result is 
that the weak Killing condition (at first order) is again sufficient to eliminate the ambiguity 
in the definition of the entropy current. In Appendix A, we continue this type of analysis 
to second order in gradients. Here one needs a new set of constraints on the components 
of the radial derivative of £ A in addition to the weak Killing condition at second order. It 
turns out that these conditions ultimately follow from the positivity of the entropy current 
divergence. 

Instead, here we propose a natural, geometrical way to eliminate all the ambiguities 
a priori and in general is to simply impose the stronger condition in ([5]). This set of 
(d + 2) equations fixes the (d + 2) components in £ A , from this one can determine the radial 
derivatives that appear in the formula ([3]). With ([5]) imposed, £ A is the vector generating 
coordinate transformations that preserve the bulk gauge condition typically employed in the 



fluid-gravity correspondence 27J] 



g rr = 0; g rfJL = Ufj, . (19) 

Geometrically, in this gauge lines of constant x M are null geodesies and r is the affine pa- 
rameter along these geodesies. This choice of gauge is associated with a trivial mapping of 
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points on the arbitrary boundary hypersurface to the event horizon along the ingoing null 
geodesies, i.e. the gauge theory coordinates are also coordinates on the event horizon 
hypersurface and r^x^). Note that the form of this gauge is the same in both the r and f 
coordinate systems, so that ultimately either choice will suffice. In the following sections, 
we will compute the perturbed Rindler solution to second order and demonstrate that the 
now unambiguous entropy current fl3]) has a non-negative divergence, consistent with the 
generalized Second Law. 

III. PERTURBED RINDLER SOLUTION IN VACUUM EINSTEIN-GAUSS- 
BONNET GRAVITY 

Following (171) the metric for the flat Rindler wedge can be expressed in the following form 

ds 2 = gABdx A dx B = —(1 + p 2 (r — r c ))u fl u u dx' 1 dx 1 ' — 2pu fl dx' 1 dr + P^dx^dx" . (20) 
The holographic fluid system lives on the surface r = r c and has a pressure 

V=^^=- (21) 

Here P^ v = rj^ + u^u V) so the fluid is on the ordinary flat Minkowski metric, > is the 
location of the horizon. Note that in this section we will use the standard radial coordinate r 
(that is in some sense adapted to the boundary) instead of the horizon adapted f. Evaluating 
the holographic (Brown- York) stress tensor at r = r c yields 

T ta ,=pP l „,, (22) 

which indicates the Rindler metric is dual to a fluid with zero equilibrium energy density 
p = 0. 

A flat metric is not only a vacuum solution to Einstein equation, but also to any standard 
higher curvature theory, where the action is constructed from higher powers of the Riemann 
tensor and its derivatives. Here we will consider the Einstein-Gauss-Bonnet theory, which 
is given by the action 

C GB = j d d+2 x^ [R + a (R 2 - AR CD R CD + R CDEF R CDEF )] , (23) 
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where a is the Gauss-Bonnet coupling constant. We consider d > 3 since for d < 3 the 
Gauss-Bonnet term is topological and does not affect the field equations. There are two 
reasons to consider the Gauss-Bonnet term. First, Einstein- Gauss-Bonnet gravity is notable 
because even though the action is higher order in the curvature, for the unique combination 
of curvature invariants in the second term of (1231) . the field equations remain second order in 
derivatives of the metric. Second, such a term often arises in the low energy limit of string 
theories. 

The field equations are given by 

Gab + 2aH AB = , (24) 

where the tensor Hab is defined as 

n-AB — titiAB — ^tiAC^B ~ lti J^ACBD + tlA J^BCDE 

- \g AB {R 2 ~ ±RcdR CD + RcdefR CDEF ) ■ (25) 



To put field equations in a more compact form, we take the trace of ( f24j) . which leads to the 
following on-shell condition: 

R-±B8. (26) 
Substituting back into the field equations, we find 

Y AB = Rab + 2aX AB = , (27) 

where 

Xab = Hab — -^abHq . (28) 

To model the Rindler fluid in local equilibrium, one allows n M (x^) and but leaves 

r = r c and the induced flat metric on it fixed. The metric is no longer a solution to the field 
equations, but one can expand and work order by order in derivatives of the fields and p. 
The metric (120]) is a solution at zeroth order, i.e. Yab = + 0(A), where A is a parameter 
that counts the derivatives of and p. The strategy for solving the equations is as follows. 
One introduces a metric correction at first order 

g = s (0) + 6gW . (29) 
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The corrected metric at first order induces a 5Y^ at the same order, which means it involves 
only radial derivatives. We want to solve for the metric Sg^ so that 

+ Y ab = , (30) 

where R^g comes from the zeroth order metric. This method can be generalized to solve 
for the metrics at higher order in A. If we have a solution to (n — 1) order g^ n ~ l \ then one 
introduces a correction 5g^ so that 

^+Y^ = 0. (31) 

Expanding out, we find 

+ 41 + 2a(6X™ + *W) = . (32) 



There is no change to the Hab tensor ( 1251) at the same order n since the curvature of the 
Rindler background is zero and any term in the variation would contain some factor of 
curvature at zero order. Thus, SX^b = 0. At first viscous order (n = 1) the results found 
previously for GR hold, because H^ B = 0- This follows just from the fact that the lowest 
order part of the Riemann tensor is first order for the Rindler solution and Hab tensor is 



quadratic the Riemann tensor and its contraction. In 25| it was shown that for any higher 



curvature theory the first order corrections vanish and the viscous hydrodynamics is the 



same as in GR. The result for the dual metric to first order is 
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ds 2 = — (1 + p 2 (r — r c ))u^Uydx tl dx 1 ' — 2pu IM dx fJ 'dr + P^dx^dx" 

+ 2p(r — r c )D(\xip)u iM u u dx^dx v — Ap(r — r^)u^P\d\ Inpdx^dx 1 ' . (33) 

One can show the general solution consistent with the boundary conditions takes the 
form 

PiWgS = -V P Ur' f P^Y^dr" (34) 



'^P a M£ = (1/2)(1 - r/r c )V^{x) - 2p [° dr' f ° dr"P x Y^ (35) 

J r J r' 

uVfyJJ = (1 - r/r c )A {n \x) + p f° dr' J° dr" (pP X<T Y^ - p~ l <S>Y r ^ - 2^^ 



(36) 
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where Y^g = R^ + 2aX A n ^. From here on we consider the second order viscous calculation, 
which is where the higher curvature corrections first appear. We only need to find Hab & t 
2nd order, which amounts to just 

rr(2) _ dWCDE A(l) 1 (0) A(l) A(l) CDEF /o 7 \ 

AB -"VA BCDE a" AB CDEF • \° ' ) 

We find 

£( r 2) = (38) 

fl£> = (39) 

3 3 
H $ = -^P 2a \ u {v d ^) u P ~ u {Ap\ u v)) + ^V 2 d^u x d\ X \u v) 

- ^p 2 pP%u^d x u u - ^p 2 d^d u u p + ^p 2 Ufl u u a x d x \np + ~P M A^ . (40) 
Since only the (fiu) components of H AB are non-zero, for the solution (|36|) . 

X<? = ( 41 ) 

X® = \u,m% (42) 

= - + ^)^ (2) g • (43) 



Apia ._„, , n _ dA ™a. „, . „„j n-L 



Hence we only need the trace # (2) g = P^ u H i2 J and P X P^H (2 J. Using the variables /C^ 



P M A P;9 (A n ff) , = P*PZd [x u a] , and = P;<9„ we find 



H [2)C c = +\p\d-2)n a ^ , (44) 



and 



^^#2 = +3p 2 fi/^A, + j^A/^ • (45) 



Putting all this together we find using f )36|) 

P^K^xl = (GRsolution) + 6ap 2 (r - r c ) ffi/fi^ + ^P M A/^ 



(46) 



u x P°5gfl = (1/2) (1 - r/r c )V; (2) (x) + (GRsolution) (47) 



u x u"6g l £ = (1 - r/r c )A( 2 )(x) + {GRsolution) + 3ap 4 (r - r c ) 2 ^ ^ O^IT^ 



(48) 



d 

To keep the equations compact we will not write the explicit form of the GR solution 
computed in 111, |l2|. The free functions A(x) and V^(x) are still arbitrary and may still 
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depend on the Gauss-Bonnet coupling constant. The functions are fixed by imposing gauge 
conditions on the stress tensor. In the Gauss-Bonnet theory, the form of the holographic 
stress tensor is modified to [28] 



= 2(K liiV - K^) + Aa{J lllv - 3 J, IU ) , (49) 

where 



36] 

= h&KK^Kl + K aX K° x K^ - 2K^ x K Xu - K 2 K^). (50) 

The first step to calculate the new stress tensor is consider the extrinsic curvature at 
second order induced by the second order metric. The result is 

^ = ^U. (51) 

As a result the pure Brown- York part of the holographic stress tensor in linear in extrinsic 
curvature has the form 

5T$ BY = l -A {2) P, u + + (GRsolution) - Upa (% x tt x „ + ^P^n^wA . 

(52) 

Next, we must consider the variation 5jjh) . From (150]) this involves the second order extrinsic 
curvature in ( 1511) times two zeroth order K^J . Due to the simple form 

rfS = ~f V* , (53) 

one can show that SJ^ vanishes identically and there is no contribution from the explicit a 
part of the holographic stress tensor. 

In the explicit a contribution to the stress-tensor, we have to calculate J)J , which involves 
two first order parts of the extrinsic curvature times one zeroth order part. Using 

K$ = ^ , (54) 

we find 

T® J = 4a (jpP^aplC* - pK^JC^ . (55) 
Putting this all together, we arrive at 

T$ = ^ m P^ + ^(A? + (GRterms) - 12pa (fi/^A, + ^ u Q a ^ 

+ 4aQpP^/C^/C a/3 -p/C MA /C^ . (56) 
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The first gauge condition we require is that 7$VP£ = 0, which fixes T/J 2) [37j . In this 
case there are no a corrections and as result (see Eqn. flUD) u x P°5gf} is the same as m 
GR. The second condition is to cancel all terms proportional to P^ u so that the pressure is 
unchanged at higher viscous orders. This fixes the in u x u u 5g^l to be 

A {2) = (GRresult) + — ^a^O^ - Ar c p 2 alC a ^ . (57) 

The remaining contributions to the stress tensor are 

T$ = (GRterms) - \2apVt x Vt Xv - AaplC^i- (58) 

Comparing to the general form of the fluid stress tensor lOj 



Tfu, = pu^Uy + pP^ - 2rjK, llu 

+ c x K\K Xv + caK^n m + c 3 tt x tt Xv + c 4 P x P°D x D a hyp 

+ c 5 JC^D\np + ceD^\npD^\np , (59) 

we see that the Gauss-Bonnet corrections only modify the c\ and C3 transport coefficients 
at second order to 

ci = — (l + 2o;p 2 ) 
p 

c 3 = ~(l + 3ap 2 ) , (60) 



with the rest of the coefficients rj, 02,04^ 



,Cq the same as in GR. These results for the 



transport coefficients are consistent with 25[, where the transport coefficients C1..4 were read 
off from the non-relativistic solution. 

IV. RINDLER FLUID ENTROPY IN EINSTEIN-GAUSS-BONNET 

We now want to calculate the entropy current in the Einstein- Gauss-Bonnet case using 
the entropy current prescription in Section II. We can express (|14p as 

4a B = -v / ^(-2AT%)- (61) 

Here we have defined 

= V M 4 - . (62) 
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Note that in the case of a Einstein-Gauss-Bonnet (and its Lovelock generalizations 29|). 
the covariant derivative of C ABCD on any index must vanish identically. In the (f, a; M ) 
coordinates 

£jr=-^+£gr , (63) 

and 

C%g = 2a (r^ - FT J" + RTt + R vf t - ^R^t\ . (64) 

Here the expressions are to be always evaluated at the horizon f = 0. 

Using the first order solution ( 1331) . we find that C^g^A = implies that £ A takes the form 
in ()6]), just an extension of the null normal to first order off the horizon surface into the 
bulk. In terms of the variables in ( 161]) we now have that 



D W = p~ 1 D\np 

) _ 



PtA® = -2Dilnp , (65) 



where the superscript in parenthesis denotes the number of gradients in the expression. 

At zeroth and first order, the Gauss-Bonnet current should also match the GR result since 
the first order solutions are the same, as we discussed above. Let's check this using ( IBTT) . At 
zeroth order, we immediately get the the GR entropy current Sq B = Attu^ 1 as expected. In 
the first order solution ( |33l) it was found previously that = and \/ — = 0. Using 

«(°) = - 
2 

Af = -pu, , (66) 
and (165]) . at first order the expression is 

4% = WCS V B m u v • (67) 

Then since 

itff = 0; R {1) = , (68) 
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on-shell, many terms in £v rur i n do not contribute at this order. The remaining term has 
the form 

= 2aRT vf{ ^ = (69) 

where we used the fact that 

with respect to components in the standard r coordinate. Thus s^q B = as expected since 
K^Uu = 0. 

Now we want to calculate the second order part of the current. It turns out that the first 
order data in ( 1651) is sufficient to elminate most the ambiguities. We just need the additional 
second order condition , which follows from (jSJ), 



*Af = 2pn^ . (71) 



Then using (IBTj) . (1631) and (1661) . the full current at second order takes the form 



+ 8irp 2 C^ f{2) u u + WnaJC^D^ In p . (72) 

The first line has no explicit a corrections, but there could in principle be Gauss-Bonnet 
corrections to the normal and metric determinant. 

First consider the metric determinant. In our standard r gauge for the metric, we have 

9tt = (73) 

9rn = ~pu^ (74) 

g^ u = P^ u + corrections . (75) 

Thus, one can show 

detg = pdet(P*P°5g Xa ) . (76) 

To second order, the transverse part of g^ v is the identity matrix plus corrections because 
P i p v9\l = ( see ©)• Using the formula det(J + A) « 1 + Tr{A) gives 

det g = p(l + P^5g$) . (77) 
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To find the correction proportional to a, consider (|46j) . In this term however the correction 

n^tlxv + ^P^apM* (78) 

is traceless, so the metric determinant to second order is the same as in GR. 
The null normal to horizon surface is defined as 

£ A = g AB d B (r-r h (x)) , (79) 

where in general the horizon location is a function of x^. At second order, the result is 

F {2) = g» r (2) | r=r _ i + ^P^d u \npD\np - ^P^dJDlnp) . (80) 
c ~? p p 

Here we have used the first order inverse metric and 

12 

r h = r c - — + — D\np+--- 81 
Thus, any a corrections will be in g^ r ( 2 \ However, this is given by 

g^) =p -lp,X u a g m ^ (g2) 

which we saw above has no corrections. Thus, we conclude that #L in Einstein-Gauss- 
Bonnet is the same as in GR. Since ^ 2 ) u ^ = ^> the first line of (J72J) is just the area entropy 
current of GR. Thus, 



u(2) _ u(2) u(2) , R o\ 
EGB ~ b GR GB ' \ 00 ) 



where (ll] |l2| 



HP 



„M(2) 
h GR 



4tt (1(k q/3 ;t*V + ^(^ a/ v) 



47T 

+ — (2Z^ InpD hip - 2D" ± D lap- P^d^K* + + ^a)£>a M • ( 84 ) 



and 



M(2) _ fi/m „2 r ^rur{2) 
GB 



8 V^gT V + 16^/0^ In p . (85) 



Thus, we just need to compute the Gauss-Bonnet parts in ( |85l) . For the details of the 
computation, see Appendix B. Ultimately we find that in terms of the fluid variables this 
has the form 

sjg = 87ra(/C a/3 /C a/3 K - 2ATTa{Q a/3 Q a ^ - l^aP^d v K vX . (86) 
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Now let's consider the divergence of the total entropy current. At lowest order the 
divergence of the entropy current is just 

= 4tt(0X)(2) = , (87) 

where in the last equality we have used the fluid equations. We can demonstrate the di- 
vergence is non-negative if we can show it still forms a perfect square at higher orders, 
i.e. 

Weqb = J + S^f = jJC"^ + 2Sjg) + 0(e 4 ) , (88) 

where S^ u is some second order tensor constructed from the fluid variables. 
At third order, the divergence takes the form 

*>A = 9^g$ + 47r(9X)(3) + d^S ~ 2KTS$ . (89) 
In the case of the diver genc e of the second order GR current, nothing changes from the 



previous calculations in 



11 



12]) and the divergence is of the perfect square form. However, 
the contribution from the divergence d^u^, which is determined by the third order fluid 
equation, u ll d v T ilv = 0, does have a corrections. In the Gauss-Bonnet case, only the 
coefficients of the fluid stress receive a corrections (16T)|) and the result is 

Att 

4tt(c>X)(3) = ~-l^ aP (( 2 - + 2K^D\np - 12ap 2 Q a x Qp X ) . (90) 

However, these a corrections are also of the required form and therefore will ultimately have 
no effect on whether the total current divergence is positive definite or not. 

Finally, let's consider the divergence of the last purely Gauss-Bonnet term above. Taking 
the divergence and imposing the ideal fluid equations yields just 

dp4$ = -16na}C^ (3n a x n px + K a x K px ) , (91) 

which again has the necessary form for positive definiteness. Putting everything together, 
we find that 



+ — (2D" ± \npD hxp - 2D\D hyp - (1 + Aap 2 )P» v d x K x + (/C£ + W x )Dl hip) 



(92) 
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and 

f-^EGB — ~p a P 

+ — (-8K, af> Dt;Dj]np + 8JC xf> Di;hxpDj hip - AD lnp/C a/ X/3 

-8/c^/c A/3 - 8^n a x n X(3 ) . (93) 

Note that the a pieces coming from the field equations exactly cancel out with pieces arising 
from the explicit a parts of the current and we are left with just the GR result. 



V. DISCUSSION 

In this work we considered the local horizon entropy current for higher- curvature gravita- 
tional theories, which in the non-stationary cases contains ambiguities related to the absence 
of a preferred timelike Killing vector. We argued that these ambiguities can be eliminated 
in general by choosing the vector that generates the subset of diffeomorphisms preserving a 
natural gauge condition on the bulk metric. While this requirement determines the entropy 
current uniquely, it is not yet clear that it guarantees in general that it has a non-negative 
divergence. As an example, we studied a dynamical, perturbed Rindler horizon in Einstein- 
Gauss-Bonnet gravity setting and computed the bulk dual solution to second order in fluid 
gradients. We showed that entropy current at second order has a manifestly non-negative 
divergence. 



Notably, our condition (jSJ) on £ A has also appeared recently in 30( outside the context 
of the fluid-gravity correspondence. Here the authors examine the approach to equilibrium 
horizon entropy using the Virasoro algebra, central charge, and the Cardy formula. One 
considers a subset of diffeomorphisms that in some sense respect the existence of a horizon 
in the bulk solution. Roughly, some of the full diffeomorphism symmetry is broken and some 
of the pure gauge degrees of freedom become physical. The authors argue that the subset 
consistent with (jSJ) is preferred and apply it to a stationary Rindler metric. They show the 
vector £ A obeys a Virasoro algebra with central extension and one can use the Cardy formula 
to count the asymptotic number of states. The result agrees with the Bekenstein-Hawking 
formula. It would be interesting to understand better the nature of (jSJ) and why it appears 
in different contexts. 
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In the fluid/gravity correspondence, condition (jSj) is necessary and sufficient to define 
an entropy current from the horizon data, which is compatible with the hydrodynamics 
dictated by the boundary stress-energy tensor. As noted above, it is an open problem to 
prove that this entropy current has in general a non-negative divergence. If true, it will 
prove in particular a Wald entropy increase theorem at least in the derivative expansion. A 
first step might be to consider the entropy current in Einstein-Gauss-Bonnet (or its Lovelock 



generalization), perhaps making use of the recent work in 3l|, |32|, which proved the Wald 
entropy increases in quasi-stationary processes. 
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Appendix A 

Here we present an alternative analysis of the ambiguities in the entropy current at first 



and second order in fluid gradients following 



24|. In particular, here one does not a priori 



impose (131) . At first order in gradients of the fluid variables, (Tl4|) contains the following set 
of unambiguous terms 

s { ut amb = % (-22™ (1) 4 0) + 4[V,£^] (1) ) , (94) 

and ambiguous terms 



s 



The "ambiguous" terms are associated with how one defines the surface gravity and radial 
derivatives of the null normal away from equilibrium. Note, that \W v O LTVr }^ = for any 
C ABCD . Using the fact that generically 

4°) = -2 A, , (97) 
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we find that this ambiguity is eliminated if 

uMf = 2k (1) . (98) 

One can show that this is equivalent to 018p at first order evaluated on the horizon. If we 
impose this condition, then to first order there is no ambiguity in the entropy current for 
any higher curvature theory of gravity. The unambiguous part of the current was computed 



in 



24J | in the case of a charged black brane background in Einstein-Maxwell theory with a 
mixed gauge-gravitational Chern-Simons term. The resulting entropy current agrees with 
the field theory calculation and has a non-negative divergence, indicating that this current 
may be the correct one. 

At second order in gradients, 

s ±1 = (-2c^a? + ^K4 2) ) - ^) • (") 

We can express £vrvr(i) j n ^ erms Q f tfl e fluid variables as follows 

£^r-(l) = a (l) u H u u + pQ.)pitu + 2U^Y$ + 7 (T^ , (100) 

where aP^ and (3^ are first order scalars constructed out of the fluid variables, Yf~ is a 
first order vector, and = P^P"d^\u a ) — \P^vd\u x . In general, these ambiguous terms 
have a somewhat complicated form relating the components of to k and it appears that 
we can say little without actually computing the perturbed solution in a particular theory 
to second order in gradients. 

In the Rindler case above, the ambiguities reduce to just 

s^l = 4ttpK^, 2) K - 8tt^ 2 V - ATipC^ uf{1) A^\ (101) 

via fl69|) . Essentially in (1 1 00 j) only the 7cr A " / term is non-zero. The first two terms proportional 
to m m above vanish if we impose the weak Killing condition at second order u^A^' = 2pK^ . 
To compute the remaining term, we can use the previous results in fl69l) and fl66|) to show 
that 

<S = SnaJC^A^. (102) 
Thus, the remaining ambiguity is in the other set of d degrees of freedom contained in 



(i) 
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If we repeat the calculation of the entropy current with fll02p . we find 

<Vgb = -167ra/C Q/3 \D^Dj hip + 3Q a x Qp X + K, a x K px - ^Afj 

- lQixa(d^ x )Di hxp - STxa(d^ x )A { l ] . (103) 

Therefore we see that the divergence will not be of the perfect square form in general 
unless we require exactly the last condition in (165|) . that P X A^ = —2D^\np. Using this 
information, we can work backwards and compute the form of £ A . The result is 

£ A dX A = 2pKfu^dx» + 2f(D^ \np)da* + df + ff (1) df, (104) 

We have parameterized the extra degree of freedom in if = ff^\ The conditions for 
positivity of the divergence of the entropy current seem to be just (d + 1) conditions on the 
(d+3) degrees of freedom in and i^ . However, we can actually calculate using the 
formula 

tv/ r = Ki r -»• -[rrv-] (1) lr-=o = (105) 

Raising (1 1041) to a vector yields 
,a d 1 



i(l + /Wf) u^ + r (ff^ + 2K^--D\np) ^ 
p ot' 1 \ p J or 



dX A 

Using this result and the first order metric, (11051) yields k,^ = p~ 1 Dlnp and thus 



(106) 



r,A 



dX A p v y cte^ dr 

Thus, up to the freedom in /W this result agrees with ([6]). Setting /^^ = is exactly what 
we need for flSJ) to hold. 

Appendix B 

In this appendix we show the detailed calculation of the Gauss-Bonnet part of the entropy 
current in ( |85l) . Using ( l64l) . we find 

L^ m u v = 2a(R t * vf ®u v - p-'RW + p~\R^u u )u^ + V 2 i? (2 V) (108) 
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Computing the Ricci tensor components and Ricci scalar is easier because we are working 
on-shell, which implies they are zero up to second order. Thus the transformation between 
the (r, x M ) and (f, x^) coordinates 

df 
dr 

J^=-^, (109) 

is trivial at second order and we can work in the original set of coordinates. The equations 
of motion (J27J) imply 

Ri 2) = 4«£( 2 )c (no) 

R ,r(2) = %fWO _ 2 ^ P ^ u °m (HI) 

dp c p Xu K ' 



The result is 



Rrt*) = zap (j-^j (a*^ Q/ V (112) 

RV) = Qap 2 ( ^-1) Q a ^. (113) 



d 

Note that these corrections already depend on a, which means they will be, in principle, of 
0(a 2 ) in the entropy current formula. Inserting these results into ( I85jl we find 

s £g) = i^p 2 aR^ f ®u v - 48a 2 p 2 f^jp) (O^^K ~ SnaK^A^. (114) 

Next, one can relate R^ rur ( 2 ) to the unbarred r components using the coordinate trans- 
formation fll09|) 

R ^ r ~ = ^w^ R ^ B > (115) 

which yields 

R^m^ = R »rvr{2) Uv _ 2p -*d X Inp^^ U u + R uX » r{1) U u ). (116) 

Using 

R r ^ x = u ^ vX + -u u W x - -u^, (117) 
(i) 2 2 y J 

we find 

R ^r{2) Uv = R^^) Uv - Zp~ 2 W x D{ lnp (118) 
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As a result, (11141) becomes 



4% = lQnap 2 R^ r ^ 



u u - 48a V 




(fi^fi^X - AtoaQ^Di lnp - MaJO^A? . 

(119) 



Finally, by direct computation using our metric solution, we find 




u u = — {-P^d p lC p x - K" X D{ lnp + 3W x Di lnp 



) 



(120) 



Putting all these components together yields (155]) and the 0(a 2 ) pieces ultimately cancel 
out. 
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